Attractive Casimir effect in an infrared modified gluon bag model 
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In this work, we are motivated by previous attempts to derive the vacuum contribution to the bag 
energy in terms of familiar Casimir energy calculations for spherical geometries. A simple infrared 
modified model is introduced which allows studying the effects of the analytic structure as well 
as the geometry in a clear manner. In this context, we show that if a class of infrared vanishing 
effective gluon propagators is considered, then the renormalized vacuum energy for a spherical bag 
is attractive, as required by the bag model to adjust hadron spectroscopy. 
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I. INTRODUCTION 

In the MIT bag model, a hadron is regarded as a finite 
region of space where quark and gluon fields are con- 
fined. The bag energy necessary to fit hadron masses 
must contain a term of the form —Z/R, where R is the 
bag radius and .Z is a phenomenological positive constant 
of order one Q . It is commonly accepted that this term 
is essentially a manifestation of zero point energies of the 
confined fields, however an understanding of this effect 
in the context of Casimir energy calculations remains in- 
triguing. 

These calculations consider gluons as free massless par- 
ticles inside the bag, an assumption based on asyrnp- 
totic freedom. Boyer, Davies, Bender and Hays 
and many others, studied massless scalar, vector and 
fermionic fields assuming they are confined in a spher- 
ical region of space. A systematic study has been de- 
veloped by Milton In all these studies, the 
renormalized zero-point energy contribution to the MIT 
bag energy contains the term but with a nega- 
tive Z (« —0.7), that is, a repulsive effect is obtained, 
instead of the attractive one necessary to adjust hadron 
spectroscopy. To the best of our knowledge, there is no 
explanation in the literature for this discrepancy. For 
a discussion of the bag model and the current status of 
Casimir calculations in the bag see refs. B and re- 
spectively. 

A resolution to the sign problem is hard to accomplish 
as in general there is no intuition whether the Casimir en- 
ergy should be positive or negative in a given situation. 
In general, the sign of the Casimir energy may depend 
on spacetime dimensionality, the type of boundary con- 
ditions, the shape of the boundaries, the spin of the field, 
etc. For a discussion of these issues, see for example refs. 



On the other hand, any treatment of QCD with bound- 
aries requires simplifying assumptions to render the cal- 
culation possible. In this regard, we note that the above 
mentioned simple model considering free massless gluons 
inside the bag only encodes nonpcrturbative information 
by imposing boundary conditions on the bag boundary, 
while the free propagator is taken over the whole range 
of momenta that can be accommodated inside the bag. 
This means, momenta ranging typically from the inverse 
bag radius up to the ultraviolet region. We would like 
to stress that infrared modes are more sensitive to the 
bag geometry, and these are precisely modes in the non- 
pcrturbative region of QCD. This situation, shows that 
any improvement to obtain an attractive Casimir energy 
must take into account intrinsic nonperturbative effects 
inside the bag such as the analytic structure QCD. 

At present, a derivation of the Casimir energy from 
first principles seems hopelessly out of reach. In this 
work we would like to propose instead a possible under- 
standing of the attractive nature of vacuum fluctuations, 
based on a simple bag model extended to incorporate an 
infrared modified gluon propagator. The associated in- 
frared behavior will be motivated by recent studies on 
the analytic structure of the confined gluon propagator 
in pure QCD O. 



II. SIMPLE INFRARED MODIFIED BAG 
MODEL 

In QCD it is believed that colored objects such as 
quarks and gluons cannot exist in asymptotic states. 
They should give place instead to the hadronic spectrum. 
These open problems are associated with a nonperturba- 
tive regime, and suppossed to be driven by the dynamical 
generation of a mass scale, which separates the pertur- 
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bative and nonperturbative phases. 

In the bag model, the boundary condition is devised 
to represent the nonperturbative physics associated with 
the finite size of colorless objects such as hadrons. 

Now, regarding the renormalized vacuum energy, if it is 
computed by considering the gluon sector only (quenched 
model) , the inclusion of additional nonperturbative infor- 
mation is a sensible modification to be considered. 

In this regard, note that such a calculation would be 
based on pure QCD proved by means of a boundary con- 
dition on a sphere of radius R. Of course, for very small 
R, the effect of the boundaries would be relevant, while, 
because of asymptotic freedom, the physics inside the bag 
would essentially be perturbative. On the other hand, for 
very large R, the model would behave as pure QCD with- 
out boundaries, displaying the intrinsic nonperturbative 
behavior of the gluon sector. 

For typical values of R associated with nonperturba- 
tive objects such as hadrons, both effects should be rel- 
evant when computing the renormalized vacuum energy: 
finite size effects representing the hadron boundaries, as 
well as intrinsic nonperturbative behavior in the gluon 
sector. In particular, possible nonperturbative analytic 
properties in this sector constitute an interesting aspect 
to be explored. 

In this section, according to the discussion above, an 
infrared modified bag model is presented. The basic 
physical inputs are the following: 

• As usual, we will consider the hadron as a finite 
region in space where the fields are confined by imposing 
a boundary condition at the bag boundary. 

• Additional intrinsic nonperturbative information will 
be parametrized in a model gluon propagator G(fc^) asso- 
ciated with an effective quadratic gluon action inside the 
bag. Although this procedure cannot be justified from 
first principles, similar ideas have been already used in 
different hadron models where some interesting results 
have been obtained This enables the computation 
of the vacuum energy by means of the usual formula. 
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lnG(fco,fc„), 



with G{k'^) in the place of the free gluon propagator 
l/(fc^ + ie). In this equation, T is the infinite time that 
the configuration exists and kf-^ refers to the eigenvalues 
of the Laplace operator in the corresponding geometrical 
configuration. As in former studies, we are considering 
a scalar Casimir energy to be multiplied by the number 
of polarization modes of the gluon field (including color) 
(see ref. 0). 

• The model gluon propagator G(fc^) incorporates non- 
perturbative effects such as pole suppression. That is, we 



will change the free 1 / (fc^ -I- ie) (fc^ = fcg — k^) propagator 
used in previous works, which is valid in the ultraviolet 
region due to asymptotic freedom, by an infrared mod- 



ified one. After a Wick rotation k 



-k\ the 



behavior of the modified propagator can be defined in 
euclidean k momentum space. 



G(-P) 



7^(P) (1/P), 

7^(0)(P)^ A>0, 



k'^ large 
k^ small. 



(2) 



This type of infrared vanishing behavior has been ob- 
tained in different scenarios (see the discussion below). 
For large fc^, we have also included a possible TZ{k'^) 
factor encoding running coupling constant information. 
Note also that the free propagator would correspond to 
A = -1 and 7^ = 1. 

When computing Casimir energies, we believe that the 
above mentioned simple model can capture the essential 
modifications implied by the analytic structure of QCD. 
With regard to the behavior in (|2J), we are motivated by 
recent progress on the form of the gluon propagator in 
pure QCD. In general, because of confinement, we know 
that the gluon propagator must suffer a dramatic change 
when A:^ — !■ 0: as gluons cannot appear in asymptotic 
states, the free pole at fc^ = must disappear in the 
complete theory. In fact, many authors coincide that the 
exact propagator should be infrared finite when k'^ — > 0. 
Indeed, an infrared vanishing behavior has been obtained 
in different scenarios: by studying the Schwinger-Dyson 
equation by restricting the path integral so as to 
avoid Gribov copies [l^ . and also in a Lattice formu- 
lation p^ . In other words, different scenarios point to 
an intrinsic analytic structure in pure QCD associated 
with an infrared suppressed propagator, describing the 
so called "confined" gluons, in contrast to an infrared 
enhanced "confining" 1/k* behavior (for a review, see 
ref. Hi). 

In particular, studies of the (Landau gauge) Schwinger- 
Dyson equation in pure QCD lead to an infrared vanish- 
ing gluon propagator G(fc^)(?7^'^ — k^k'^/k'^), with the 
form 

G(P) = (-p _ ie)\-k^ + A2 - ^e)-(^+l)7^(-P - ze), 

(3) 

where A > and A is of the order of Aqcd- Note that 
after a Wick rotation, —k^ — it fc^, this ansatz sat- 
isfies the infrared vanishing behavior given in (O. In 
the Schwinger-Dyson context, the additional condition 
A < 0.4 is satisfied, and the factor TZlk^) represents the 
perturbative running. That is, for k^ >> A^, 



7^(P) 



const. 
ln(P/A2) 



7 > 0. 



(4) 

On the other hand, the Gribov mechanism is based on 
the restriction of the path integral domain so as to avoid 
copies of the gauge fields in the pure gluon theory defined 
on euclidean space. This is necessary as in the nonper- 
turbative region the usual gauge conditions do not fix 
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the gauge completely. In Landau gauge, the associated 
nonperturbative gluon propagator is given by, 



GGribov(-^^) 



/c4 + M4 ' 



(5) 



which satisfies the behavior (0) with A = 1 and TZ = 1. 

As we are concerned with a model gluon propagator, 
we will first develop a Casimir energy calculation based 
on the associated general analytic properties displayed 
in Minkowski space, not relying on a specific form for 
the propagator. We will obtain a representation for the 
Casimir energy where the analytic structure of the model 
and the effects of geometry will be clearly distinguished. 
Only at the end we will discuss the effect of pole suppres- 
sion on the renormalized vacuum energy. For this aim, 
we will consider for instance a model gluon propagator of 
the form given in eq. which displays a well defined 
analytic structure in Minkowski space. In this regard, 
note that the Gribov ansatz cannot be rotated back to 
Minkowski space because of the complex poles in eq. 
as is well known, the Gribov scenario is only defined in 
euclidean space. 



III. CASIMIR ENERGY AND THE ANALYTIC 
STRUCTURE OF G 

In order to study the effect implied by the analytic 
structure of the modified propagator, we will start by 
representing the logarithmic derivative of G in terms of 
the decomposition, 

d f°° 1 

^lnG(^ + .)^y^ ^^W)-^^:^, (6) 

where A= and /3 is given by, 

(7) 

Defining the jump at the discontinuity of an analytic 
function F{A), 



S[F{A)] = F{A + ie) - F{A - ie), 



we can also write, 



-InGiA) 



(8) 



(9) 



In this manner, from eq. ((SJ, we have the representation, 

\nG{A + ie)^C+ / d^^/SOi^) In G^2 + ie), (10) 
Ja 

with C — const, and G^2(.A) = {—A + /i^)~^ being the 
free propagator for a field with mass parameter /it^. 

In fact, when defining the Casimir energy in eq. a 
parameter with the dimension of [mass]^ multypling 



G must be introduced, so as to have a dimensionless ar- 
gument in the logarithm. In an equivalent manner, we 
can introduce a parameter [^] = mass, multiplying 
G^2 in the second member of eq. UlUI) . and absorb the 
constant C in its definition; then, we have, 



/•OC 

In(x'G) = / V/3(A*')ln(^'G, 
Ja 



(11) 



Now, considering this replacement in eq. 1^, we can 
take the formal trace in the second member of eq. Hll|l . 
to obtain the Casimir energy in terms of the regularized 
representation. 



2tt 



In other words we can write. 

En = / dt?P[t?)Ec{li''). 



rog 



where 

i?c(/i') 



E 



"2^ 



In {eG^2{An+i€)) 



(12) 



(13) 



(14) 



is the Casimir energy for free field modes with mass pa- 
rameter /i^, obtained by means of the regularized par- 
tition function {An — kg — k'^). After a Wick rotation 
ko — > iko, we can also write. 



kg + fc^ 



(15) 



rog 



which can be defined, for instance, by using the zeta func- 
tion rcgularization technique. 



^ I dko fk^ + kl+ fi^"- 



(16) 



(17) 



As discussed in ref. [I^, the zeta function result for 
Eci^J^) in eq. H16|l differs in general from the Casimir 
energy obtained by the zeta function rcgularization of 
the summation of the free field energy modes, -Emodc = 

(l/2)E„(fc^+A^')'/'ke„ 



i?modc(M') = i lim \ [C3(-l/2 + e) + C3(-l/2 - e)] 



C3(s) 



kl + 



(18) 
(19) 
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In that reference, the following relationship has been es- 
tablished, 



Ec{^^^) = i?modc(M') + I (V'(l) - ^(-1/2)) 



47r2 



(20) 

where ip{s) is the digamma function and C2 is the second 
Seeley-De Witt coefficient. 

While the C2 coefficient for a geometry consisting of 
two parallel plates vanishes, in the case of a spherical bag 
it is nonzero. As noted in ref. this difference reflects 
the inherent ambiguity introduced in the Casimir energy 
when removing the polar part in cq. (|f 9|l . by choosing 
the principal value, such as in eq. H18|) . or by using other 
possible prescriptions. 

According to ref. C2 can be read from the varia- 
tion of -Emodo under a change of the parameter ^, 



mode 1^' 



mode 1^ 



47r2 



(21) 



We will discuss below the particular form of C2 in the 
bag model context, and we will see that both Casimir 
energy definitions, Ec and -Emodc, lead in fact to the 
same physical implications. 

It is important to remark that in the superposition we 
have derived for the Casimir energy Ec (cf. eq. H13|l ). 
the Casimir energies Ec{y''^) encode the geometry, while 
the factor encodes the analytic structure of the 

model. Note also that the propagator for a free massive 



field of mass m would give (3{fi ) = d{p, 



-m 



in eq. lO, 



and eq. ((T^ would give Ec — Ec{m^), as expected. 

In general, because of the presence of Ecifi'^), our ex- 
pression corresponds to a physical situation where all the 
/i^ field modes that participate in the representation © 
are confined, that is, they see the boundary condition. 
From this point of view our approach is quite natural, 
since confinement is supposed to act over all the gluon 
modes represented by the model propagator G. This 
should be contrasted with typical radiative QED correc- 
tions to the Casimir effect where, in the computation 
strategy, the boundary condition is seen by the free pho- 
ton propagator but not by the continuum of electron- 
positron modes represented in the vacuum polarization 
[l9|. In these cases the overall result is a renormaliza- 
tion of the sphere radius, the distance between plates, 
etc. This renormalization can be traced back to the ef- 
fect of electron-positron pairs, which are not confined by 
the boundary conditions, thus leading to an effective en- 
largement of the confining region. 

We also note that since our formula relies on the su- 
perposition of the usual Casimir energies, obtained from 
the regularized partition function for field modes with 
mass parameter /i^, the form of the divergences we will 
obtain here are the standard ones. For a detailed dis- 
cussion see refs. [23,0,122. These divergent terms will 
renormalize similar terms in the bag model to be fixed 
by phenomenology. For instance, the divergent volume 
term will renormalize the bag constant B. 



Now, according to our discussion at the end of section 
iim let us consider for the model gluon propagator an 
ansatz of the form given in eq. Q. In this case, eq. © 
leads to. 



/3(fc') 



1 

27ri' 



A 



1 + A d\nn{A) 



fc2- A2 



dA 



^XSiP) + (1 + X)6iP - A2) + (3niP), 



(22) 



where we have used eq. (O, the property 
and the definition, 



1 

27ri' 



d\iin{A) 



dA 



(23) 



Then, from eqs. and 1221) . the Casimir energy for 

the infrared modified model results. 



Ec - -XEc{0) + {^+X)Ec{A')+ / dti'f5u{^^')Eciti'). 

Jo 

(24) 

Of course, if a free massless propagator were considered, 
that is, A = —1 and TZ = 1 (cf. eq. lO), we would have 
/3(/x2) = S{^'^), and the Casimir energy Ec — Ec{0) for 
a massless field would be reobtained. 

In that case, the zeta function calculation for the reg- 
ularized energy mode summation in eq. (|I8() has been 
performed in ref. pQ]. obtaining. 



Err 



1 

R 



0.08392 



3157r 



(25) 



which corresponds to photon modes confined inside the 
bag. Using this expression in eq. (|^ . we read. 



47^2 



1 



3I57ri?' 



(26) 



which corresponds to an additional 1/R term when pass- 
ing from the Casimir energy defined by i?modc(0) to 
Ec{0) in eq. if^ : evaluating the digamma functions, 
we get. 



Ec{0) 



1 

R 



0.08640 + 



3157r 



(27) 



In ref. , the zeta function result for the summation of 
energy modes in eq. I|25|) has been compared with that 
obtained in ref. Iji 



1 

R 



0.08984-1- 



SIStt 



ln(5/8) 



(28) 



found by the Green's function method, where (5 is a cut- 
off associated with a nonzero "skin depth" representing a 
realistic boundary, instead of a sharp mathematical one. 
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As noted in that reference, the 1/R parts obtained from 
the first term in eqs. (|25|l and (|28|l . differ within a 6.7%, 
and in fact, there is no reason why these parts should be 
equal, as they may vary by just changing the values of 
the parameters 5 and ^. The same situation applies to 
the 1/R part coming from Ec^ obtained from the first 
term in eq. (|27|l . 

In a bag model context, the point is that a natural 
scale exists which permits to derive some qualitative and 
semiquantitative conclusions (see ref. and references 
therein). 

For instance, the plausible value for ln(5/8) is of or- 
der one, when realistic boundary conditions with a "skin 
depth" of about 10% of a typical bag radius is consid- 
ered, see ref. There, a first crude estimate ~ /R 
has been obtained for the zero point energy by dropping 
the logarithm term in eq. H28(l (which corresponds to 
10% of the first term), and multiplying by a factor eight 
counting for the number of gluons. As stated in 9], it is 
certainly very hard to doubt about the sign of the effect. 

In the zeta function regularization scenario, according 
to refs. |23|, and 0|, in a pure QCD context ^ should 
be associated with the energy scale parameter Aqcd', as a 
consequence, we will see that the discussion in the above 
paragraph is also realized in this type of scenario. 

Let us consider ^ = A = 600 Mev (a typical value for 
the fitting of the nonperturbative propagator Q), and 
the radius R taking values around Rq = Mev^^ (a 
typical value for a hadron radius). We see that for this 
choice the logarithm in eqs. H25() and (|27(l is of order 
one, and again the corresponding second term represents 
about 10% of the first term. In general, we can write eq. 
(P7|l as. 



Ec{0) 



1 
R 



0.08640 ■ 



3157r 



ln(Ai?o) 



3157r 



(29) 

After multiplying by eight colors, the first two terms in 
the bracket give a contribution +0.75/i?. Then, we see 
that in the zeta function regularization scenario, a clear 
scale for the Casimir effect again arises. The Casimir 
energy in the "photon-like" bag context can be consid- 
ered as consisting of a ^ +0.7 /R part, plus logarithmic 
corrections. 

We also note that even in an R interval that covers two 
orders of magnitude around Rq, from (1/10) Rq to 10 Rq, 
the whole bracket in eq. (|29|l remains of order one, vary- 
ing from +0.6 to +0.9. Then, in this scenario it is also 
very hard to doubt about the repulsive character of the 
effect. For "photon- like" gluon modes confined inside the 
bag, in order to change the sign of the Casimir energy and 
comply with phenomenology, the typical hadron radius 
Rq and typical scale A in pure QCD should lead to a neg- 
ative contribution of order one coming from the first two 
terms in eq. H29|l . This would happen for ARq < 10~^, 
which is an unrealistic situation. 

Coming back to the infrared modified bag model, we 
observe that the first term in eq. (|24|l only involves the 



soft /i^ = modes. On the other hand, note that the sec- 
ond term in eq. (|24|l is associated with a massive Casimir 
effect Ec{A'^) which is expected to be suppressed. 

This can be seen by considering the general form of the 
Casimir effect for a massive scalar field. In ref. 22|, the 
zeta function regularized expression for the Casimir en- 
ergy has been computed; the associated divergent terms 
and ambiguities always display nonnegative powers of the 
mass of the field. On the other hand, as discussed in ref. 
|2lj . in the massive case there is a natural renormaliza- 
tion prescription, requiring that in the infinite mass limit 
all effects coming from quantum fluctuations should van- 
ish. In this manner, after renormalization, it is obtained. 



fjAR) 
R 



(30) 



where / is a well defined, finite and unambiguous func- 
tion. 

Although a simple analytic expression for f{AR) is 
lacking, its numerical evaluation has been performed in 
ref. [22|, and the result for the interior problem with 
Dirichlet boundary conditions has been presented by 
plotting / = REc as a function of the dimensionless 
variable AR. By using this information, it is simple to 
estimate the second term in eq. (|24|l . taking as before 
A — 600 Mev, and R = or R = ^ (typical val- 

ues in Mev~^ for the pion and proton radius, respec- 
tively). Around these values, we have Ec{A'^) = f{2)/R 
or f{3)/R, respectively. From ref. |22l |. we can ex- 
tract the estimate /(2) w .0002 which, after including 
a 2 X 8 factor counting the number of physical polariza- 
tions times the number of colors, implies that the second 
term in eq. (|24|l typically corresponds to a 1% of the first 
term (A is of order 1). The analysis for /(3) gives an even 
higher suppression than the previous one. 

Note that in our model the first two factors in eq. Q 
are the relevant ones to interpolate between an ultraviolet 
asymptotically free and an infrared vanishing behavior of 
the effective gluon propagator. So that, in principle, the 
model could be defined with TZ = 1 and the third term 
in eq. H24|l would be absent. 

However, it is interesting to analyze the effect of a non- 
trivial TZ factor. To be more specific, let us consider for 
instance the ansatz for TZ — a~'^ (cf. eq. 01), now ex- 
tended over the whole range of momenta according to 



TZ^ i^A- A^^o!-^ , 



A-TT 

-a(0)A2 + —A 



A2 



ln(-^/A2 



A + A^ 



(31) 



, (32) 



where a(0) = 8.915/iVc, 7 = {-UNc + ANf)/{22N^ - 
ANf), b = (lliVc - 2iV/)/3, Nc and Nf being the number 
of colors and fiavors, respectively, and the ie prescrip- 
tion in the variable ^ = A;^ is understood. This form 
implies a cut in Minkowski space along the timelike fc^ 
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axis with branch point at fc^ = 0, and the contribution to 
the discontinuity (3n{^^) defined in eq. (|23|l has support 
starting at /z^ = 0. 

Using eq. H23I) . we have, 



d_. 
dA 



(33) 



The first term in f3fi simply renormahzes the second term 
in eq. (|23, that is, 



Ec - -XEc{0)+ 

+ (l + A-7)^c(A') + 



(34) 



and as 7 is of order one, the same analysis we have done 
before can be applied to conclude that the second term 
in eq. is suppressed. 

On the other hand, can be written in terms of the 
real and imaginary parts of a' {A + is) , 



Pni^'^) = ^Tjarctan 



TT dA 



3(a') 



(35) 



A-TT 

5R(a') = -aiO)A^ + —A 



ln(^/A2 



A2 



ln^(^/A2) +7r2 A + A^^ 



9(a') = 



A 



b In2(^/A2) + ^ 



(36) 



(37) 



In Fig. n we plot PjiifJ.^) as a function of fi. We observe 
that the function displays a definite sign, weighting 
the massive Casimir energies in the third term of eq. (|34|l . 

The total weight that is distributed over the whole 
range of masses is, 



7 /5R(a') 
- arctan — — 



-7, (38) 



which can also be verified numerically; for Nc = 3, and 
Nf = (quenched theory), we have —7 « 0.59. 

Now, we can use the general parametrization of the 
Casimir energy in the massive case, given in cq. (|30|l . 
with A — > /i, to write, 



dfi'(3'T^{fi')Ec{ti' 



d^l'^3'^^{^l')f{^lR). 



(39) 

In order to obtain a bound for the contribution of this 
term, we can consider, as before, a typical radius R — 
2^ Mev~^, and divide the integration in three inter- 
vals i = Q,l, 2: from to (200 Mev)^, (200 Mev)^ to 
(400 Mev)2 and (400 Mev)^ to 00. 

Note also in Fig. ^that in the limit /i^ 0, we have 
(^n ~^ 0- Indeed, it can be verified that tends to zero 




FIG. 1: This is a plot of (3^ (in Mev~^), given in eq. I|85^ . vs. 
the mass parameter (in Mev). We have considered Nc — 3, 
and Nf = (quenched case). 



like ~ 1/ (In ^2)2^ Therefore, in eq. (|39|l . the contribution 
coming from a region very close to fi^ — is suppressed. 
Outside from this small region, it has been shown in ref. 
[22I that f{fJ,R) is positive definite, displaying a maxi- 
mum value /o ~ .0030 in the region Iq. In the region 
/ii? > 1, the function f{fJ,R) decreases monotonically to 
zero. Then, the maximum values in regions Ii and I2 
are given by /i = /(I) « .0005 and = /(2) « .0002, 
respectively. These values have been estimated from the 
numerical result presented in that reference. 

Then, replacing in each interval /(/ii?) by fi, we ob- 
tain, 



/"OO 2 „ 

^0 Jh 



(m^) (40) 



With regard to the integrals in eq. 1001), we obtained 
.001 and .006 in the intervals Jq and Ji, while in the 
interval I2 the integral completes the total weight .59. In 
this manner, around a typical radius of R 



200 



Mev" 



d^l'(3'^^{^l')Ec{^i') < io^Vr. 



(41) 



After including a factor for the total number of polar- 
izations, we see that this term represents a contribution 
less than 0.3% when compared with the first term in eq. 

As expected, the effect of the boundaries on the finite 
part of the zero point energy is dominated by the soft 
modes iJ? = present in the first term of eq. 1)34(1 . We 
have also considered a number of flavors Nf — 3, obtain- 
ing no essential modifications in the analysis. 

Then, in the context of the infrared modified bag 
model, the 1 /R part of the renormalized vacuum energy 
turns out to be —Z/R, with Z ^ -1-0. 7A, which is attrac- 
tive for an infrared vanishing model propagator (A > 0). 
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IV. CONCLUSIONS 

In this article, we have computed the Casimir energy 
in a bag model containing a modified gluon propagator. 
Our motivation comes from previous attempts to describe 
the vacuum energy in the bag in terms of Casimir energy 
calculations. These attempts, based in a model where 
the gluons are associated with massless fields confined 
on a spherical region, failed to describe the attractive 
nature of the vacuum energy necessary to adjust hadron 
spectroscopy. 

In the modified framework we have considered, we were 
able to clearly trace the effects introduced by the analytic 
structure of the model, on the one hand, and the geome- 
try of the boundaries, on the other (see equation ifT^ '). In 
this context, we can see that the introduction of pertur- 
bative information in the above mentioned massless gluon 
models is not enough to render the finite part of the zero 
point energy attractive; if only a running coupling were 
taken into account, the effect would be suppressed when 
compared with the main Casimir effect coming from the 
pole at fc^ = 0, which is repulsive. 

Then, in order to understand the attractive nature of 
the Casimir energy due to gluons, the consideration of in- 
trinsic nonperturbative effects is fundamental, however, 
an analytic approach from first principles is a formidable 
challenge. For this reason, we have introduced the sim- 
ple infrared modified bag model with an effective gluon 
propagator. Within this context, we have seen that the 
repulsive or attractive nature of the vacuum energy de- 



pends on whether the model propagator describes "con- 
fined" gluons with infrared vanishing behavior, like in the 
Schwinger-Dyson and Gribov scenarios, or it is singular 
at = 0. This last situation is verified, for instance, 
when a free or a "confining" l/k'^ model gluon prop- 
agator is used. 

From this point of view, we see that an infrared vanish- 
ing gluon propagator is preferred as it corresponds to an 
attractive term ~ — 0.7A/i?, A > 0, in the renormalized 
vacuum energy. Hadron phenomenology requires a gluon 
contribution of the form —Z/R, with Z of order one (see 
the discussion in ref. p^l. In the simple model we have 
presented, this would be achieved with an infrared behav- 
ior of the model gluon propagator of the form ~ (A:^)^, 
with A of order one, which corresponds to a typical or- 
der of magnitude appearing in those pure QCD scenarios 
where an infrared vanishing propagator for "confined" 
gluons is obtained. 

As further studies suggested by the present work, we 
point out the consideration of similar nonperturbative ef- 
fects in the fermion sector, as well as in other bag models 
where the hyperfine hadron structure is analyzed |24| . 
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